Atmospheric neutrino oscillations and tau neutrinos in ice 
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The main goal of the IceCube Deep Core Array is to search for neutrinos of astrophysical origins. 
Atmospheric neutrinos are commonly considered as a background for these searches. We show here 
that cascade measurements in the Ice Cube Deep Core Array can provide strong evidence for tau 
neutrino appearance in atmospheric neutrino oscillations. A careful study of these tau neutrinos is 
crucial, since they constitute an irreducible background for astrophysical neutrino detection. 
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I. INTRODUCTION 

Over the last decade, a large number of experiments of 
different types have provided strong evidence for neutrino 
oscillations and thus for physics beyond the Standard 
Model, see Ref. [l| and references therein. 

Cosmic ray interactions in the atmosphere give a nat- 
ural beam of neutrinos. These atmospheric neutrinos in 
the GeV range have been used by the Super-Kamiokande 
(SK) detector to provide evidence for neutrino oscilla- 
tions [2| . The large size of neutrino telescopes such as 
AMANDA, IceCube and KM3NeT makes possible the de- 
tection of a large number of atmospheric neutrino events 
with a higher energy threshold, ~ 100 GeV, even though 
the neutrino flux decreases rapidly with energy (~ E~ 3 ). 
Built to detect neutrinos from astrophysical sources, or 
from the decay of Weakly Interacting Massive Particles 
(WIMPs) annihilations [3[ , at the high detection thresh- 
old energies of these ice/ water Cherenkov detectors, neu- 
trino oscillation effects would be small. 

Recently, a low energy extension of the IceCube de- 
tector, the Ice Cube Deep Core array (ICDC) has been 
proposed and deployed [J] . It consists of six densely in- 
strumented strings (7 m spacing among optical modules) 
located in the deep center region of the IceCube detec- 
tor plus the seven nearest standard IceCube strings. Its 
goal is to significantly improve the atmospheric muon 
rejection and to extend the IceCube neutrino detection 
capabilities in the low energy domain, down to muon or 
cascade energies as low as 5 GeV. The instrumented vol- 
ume is 15 Mton. Such a low threshold array buried deep 
inside IceCube will open up a new energy window on 
the universe. It will search for neutrinos from sources in 
the Southern hemisphere, in particular, from the galac- 
tic center region, as well as for neutrinos from WIMP 
annihilation, as originally motivated. In [5J we have pro- 
posed neutrino oscillation physics as a further motiva- 
tion for building such an array. In particular, we have 
analyzed the sensitivity of ICDC to the neutrino mass 
hierarchy. ICDC can detect up to 100, 000 atmospheric 
neutrino events per year, orders of magnitude beyond the 
present data sample, providing rich opportunities for de- 
tailed oscillation studies. In the same spirit of Ref. @, 
we concentrate here on the neutrino oscillation analysis in 



ICDC, focusing now on the cascade signal. By exploiting 
the cascade channel, ICDC could firstly provide strong 
evidence for tau neutrino appearance from oscillations 
of atmospheric neutrinos, greatly improving previous SK 
results on tau neutrino appearance evidence @. 

Section [TT] reviews briefly our present knowledge of the 
neutrino oscillation parameters as well as their expected 
errors from near future facilities. We then describe the 
cascade analysis details, presenting our main results in 
section IIIII 



II. NEUTRINO OSCILLATIONS AND 
CASCADES IN ICDC 

Neutrino data from solar, atmospheric, reactor and 
accelerator experiments is well understood in terms of 
three-flavor neutrino oscillations. Two Am 2 values and 
two (large) mixing angles are well determined, while the 
third mixing angle is limited to be very small. The CP- 
violating phase (8) is completely unconstrained. In addi- 
tion, the sign of Am| x is also unknown. 

The best fit oscillation parameter values obtained from 
present data are Q: 
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and sin 2 26»i 3 < 0.15 for Am^ = 2.5 x 10~ 3 eV 2 . No- 
tice that an extra unknown in the neutrino oscillation 
scenario is the octant in which 6*23 lies, if sin 2 2023 1- 
This has been dubbed in the literature as the 023 octant 
ambiguity. 

In the near future, long baseline experiments like MI- 
NOS [8| will improve the current precision on Am 2 ^ 
and possibly discover a non-zero value of 0i3, if this is 
close to the present upper limit. In a few years, reac- 
tor experiments like DoubleChooz Q, RENO M and 
DayaBay [ll| will provide improved sensitivity to 0i3. 
This information can be used as input in our analysis, 
reducing some of the parameter uncertainties. 
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FIG. 1: survival probability and f M —¥ v T oscillation proba- 
bility as a function of the neutrino energy (in GeV), assuming 
upward going neutrinos and sin 2 2613 = 0.1. 

In the past, atmospheric neutrinos in the SK detector 
have provided evidence for neutrino oscillations and the 
first measurements of lAmfjJ and sin 2 2^23 0) therefore 
providing compelling evidence for neutrino oscillations 
versus other more exotic phenomena 12]. While facing 
more systematics than accelerator/reactor experiments 
due to the uncertainties in the (natural source) neutrino 
fluxes, atmospheric neutrinos provide great opportunities 
for exploring oscillation physics due to the large range of 
energies and pathlengths that they span. The ICDC de- 
tector will collect a data sample which is a few orders of 
magnitude larger than the SK experiment and can also 
measure energy and directional information, such that 
many of the systematic errors associated with unknown 
normalizations of fluxes, cross-sections, etc. can be much 
better understood by using the data itself. 
The IceCube detector and its Deep Core extension are 
optimized for detecting muon tracks from the charged 
current interactions of v^. It is, however, possible to 
also detect cascades |13j . While these type of events no 
longer provide directional information, their energy can 
be measured quite precisely. We analyze here the cascade 
events in the ICDC array. There are several contributions 
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where tit is the number density of targets, V is the vol- 
ume of the detector, 6 is the zenith angle direction of the 
neutrino and t is the observation time. The detector is 
a cylinder of 250 m diameter and 350 m height, i.e. the 
total physical mass is of around 15 Mton. The low en- 
ergy events however are mostly single-string events and 



to the cascade signal: charged current interactions of v e , 
neutral current interactions of all neutrino flavors and 
electromagnetic and hadronic decays of tau leptons pro- 
duced in charged current interactions of v r . There are 
several important observations which suggest that the v T 
signal can be significant and can provide evidence for v r 
appearance from oscillations of v^: 

• First, the atmospheric electron neutrino and electron 
antineutrino fluxes at the relevant energies are signifi- 
cantly lower than the muon neutrino flux, such that the 
v e charged current interactions do not completely over- 
whelm the event rate. 

• In addition, the energy range covered by ICDC corre- 
sponds to a maximum of — > v r oscillations (minimum 
of survival), as can be noticed from Fig. [TJ 

The large flux of atmospheric muon neutrinos can thus 
lead to a large flux of tau neutrinos. This has already 
been noted in where it was shown that — > v T — > 
t — y f± provides a non-negligible contribution to the muon 
track rate. It is also important to note that, unlike for the 
SK detector, which is sensitive at much lower energies, 
tau threshold production effects are relatively small, only 
affecting the lowest energy events detected by ICDC. 

We investigate the neutrino energy range between 
10 GeV and 100 GeV, assuming bins of 5 GeV width 
in the observable energy. Cascades have very little di- 
rectional information, especially at these low e nerg ies, so 
we integrate over all upward going directions [19]. The 
downward going neutrinos are largely unaffected by os- 
cillations, so they can be used for determining the atmo- 
spheric neutrino flux and thus the contribution of the v e 
charged current interactions to the overall cascade rate. 
In our numerical calculations we have taken into account 
full three flavor oscillations. It is however straightforward 
to see that solar parameters do not play an important role 
in the analysis due to the rather high energy threshold 
of ICDC. Also, 6*i3 effects, while in principle observable 
for values of #13 close to the present bound, do not affect 
any conclusions regarding v r rates, which are determined 
by the (maximal) atmospheric mixing angle #23- 

The spectrum of v e induced events as a function of 
observable energy is given by: 



■XEu, o) + d ^ d n ] p ^ +("-►*). (5) 
I 

we consider the effective volume in this case to be six 
cylinders of 40 m radius centered around the six densely 
instrumented ICDC strings. 

The second term of Eq. ([S]), which contains the con- 
tribution from oscillations of — ¥ v e , is negligible in 
practice. This contribution is very small even for the 
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maximum current allowed value of sin 2#i3 and after 
considering enhancements in the oscillation probabilities 
due to matter effects inside the Earth. It could, however, 
become relevant in the presence of non-standard neutrino 
interactions. 

For charged current electron neutrino interactions, the 
entire energy of the incident neutrinos is transferred to 
the cascade and thus measured. 



For neutral current interactions, only a fraction y of the 
initial neutrino interaction is transferred to the cascade. 
Due to the steep energy dependence of the atmospheric 
neutrino fluxes, their contribution is thus expected to be 
smaller. The spectrum for neutral current (NC) inter- 
actions as a function of the observable energy is given 
by: 
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For tau neutrinos, charged current (CC) interactions tromagnetic and hadronic, depending on the tau decay 
lead to two types of contributions to cascade events: elec- mode. The rate for hadronic events is given by: 
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For our numerical calculations, we use differential neu- 
trino cross-sections as given by [lij with CTEQ6 par- 
ton distribution functions. In terms of dimensionless 
variables typically used for these cross-sections we have: 
da^ c /dEr ~ l/E v da^ G /dy with y = 1 - E T /E U and 
dn had /dE T ~ E 2 ons /E T dn had /dz with z = l- E obs /E T . 
It is important to note that for tau neutrinos, tau thresh- 
old suppression is still important in the lower energy 
range discussed here. We have included the threshold 
corrections following Ref. [Hj]. The 3.5 GeV tau lepton 
production threshold energy has made the tau neutrino 
detection very difficult up to now, since the atmospheric 
neutrino flux for energies above tau lepton production 
threshold is very low for existing detectors. For instance, 
for the SK experiment, and assuming maximal mixing in 
the — v r sector, only one CC tau neutrino event is 
expected per kton-year of exposure. The ICDC exper- 
iment, which benefits from a much larger instrumented 
volumen than the SK experiment, is in an unique posi- 
tion to detect tau neutrino interactions, as it has good 
sensitivity in an energy range where neutrino oscillations 
lead to a large number of tau neutrinos and, at the same 
time, is high enough to no longer be strongly affected by 
tau threshold suppression. 



The rate for electromagnetic cascade events is given by 
an expression similar to Eq. ((7|) , the only difference being 
in the decay rate dn/dE T . 
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FIG. 2: Electron CC and all flavors NC events for one year 
ICDC exposure, see text for details in the calculation. 



III. RESULTS 

Figure [5] shows the cascade rates coming from v e 
charged current interactions and neutral current inter- 
actions of v e and v^. These events are indistinguish- 
able from the v T cascade events that are of interest to 
us, so they constitute an irreducible background for the 
v T search. Figure [3] illustrates the cascade rates from 
charged current v T interactions followed by hadronic or 
electromagnetic tau decays, as well as v T neutral current 
interactions. 

While the background is significant, the number of 
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FIG. 3: Tau cascade events for one year of ICDC exposure. 



events is very high. Neglecting systematics errors, the 
statistical significance could be defined as 



S = 



y/n s + n b 



(8) 



n s being the number of v r events and rib the number 
of background events (from v e changed current interac- 
tions and v e41 neutral current interactions). A statisti- 
cally significant (3a) u T appearance signal could be ob- 
tained in only a few months of observation. However, 
systematic uncertainties will limit the analysis. One of 
the biggest uncertainties is the energy threshold for cas- 
cade identification. At the lowest energies considered 
here, muon track events would be very short [2(| and 
therefore, cascade events and muon track events would 
be indistinguishable. If cascade events can not be dis- 
tinguished from muon tracks, the background to the v T 
signal becomes significantly higher (by about an order 
of magnitude) due to the contribution from the muon 
event tracks. Even in the case where muon track events 
and cascade events can not be distinguished at low en- 
ergies, the v T signal can become statistically significant 
in a year of exposure. Good muon track reconstruction 
would be extremely useful for reducing this systematic 
uncertainty and might be achievable with the planned 
additional strings in the center of the detector. There 
are other systematic uncertainties affecting the analysis, 
as the knowledge of the interaction cross-sections, atmo- 
spheric neutrino fluxes, and other neutrino oscillation pa- 
rameters. However, these systematic errors are expected 
to be under control in the next few years, exploiting data 
from future reactor and accelerator experiments as well 
as atmospheric neutrino data from the ICDC experiment 
in different angular and energy ranges than the ones used 
for the v T analysis. Appearance of tau neutrinos from os- 
cillations of atmospheric is thus likely to be detected 
in the near future. 

IV. OUTLOOK 

The IceCube detector and its Deep Core array provide 
a great opportunity for studies of atmospheric neutrinos. 



Being the largest existing neutrino detector, it will ac- 
cumulate a huge number of atmospheric neutrino events 
over an enormous energy range, thus allowing for detailed 
studies of oscillation physics, Earth density, atmospheric 
neutrino fluxes and new physics 16]. In order to extract 
all this information it is necessary to use energy and an- 
gular distribution information, as well as flavor compo- 
sition, all possible to obtain with the IceCube detector. 
Qualitatively, there are three main energy intervals and 
three main angular regions which are sensitive to differ- 
ent types of physics. 

At very high energies, above 10 TeV, neutrino inter- 
action cross-sections become high enough that neutrinos 
going through the Earth start getting attenuated [l7j . 
This effect is sensitive to neutrino interaction cross- 
sections and to the density profile of the Earth. 

The "intermediate" energy region, between 50 GeV 
and 1 TeV can provide good information about the at- 
mospheric neutrino flux, which can be used to improve 
the uncertainties in the simulated atmospheric neutrino 
fluxes [3. 

In our paper we concentrated on the "low" energy re- 
gion, below about 40 GeV, where neutrino oscillation ef- 
fects can be significant. Although the IceCube detector 
has higher energy thresholds, its low energy extension, 
the Ice Cube Deep Core array (ICDC), extends the Ice- 
Cube neutrino detection capabilities in the low energy do- 
main, down to muon or cascade energies as low as 5 GeV. 
We have shown that tau neutrinos from oscillations of at- 
mospheric muon neutrinos could be detected in the next 
few years with a high significance level, providing there- 
fore the first sizable sample of tau neutrinos. This large 
number of tau neutrinos could allow for measurements of 
v T interaction cross-sections and studies of non-standard 
neutrino interactions, which are largely unconstrained in 
the tau sector. In summary, ICDC offers a unique win- 
dow towards a better understanding of neutrino proper- 
ties due to its very high atmospheric neutrino statistics. 
Careful studies of the expected atmospheric neutrino os- 
cillation signals in ICDC, as the one carried out here 
for the v T appearance signal, are mandatory, since atmo- 
spheric neutrinos constitute an irreducible background to 
astrophysical neutrino searches. 
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